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, We comment on the reweighting method in the chemical potential (fiq) direction. 

f— *^ ■ We study the fluctuation of the reweighting factor during Monte-Carlo steps. We 

^v^j ' find that it is the absolute value of the reweighting factor that mainly contributes 

to the shift of the phase transition line (/3 C ) by the presence of fi q . The phase 
^s^j , fluctuation is a cause of the sign problem, but the effect on /3 C seems to be small. We 

i also discuss /3 C for Iso-vector chemical potential and /3 C determined from simulations 

with imaginary chemical potential. 

Qh, 1- Introduction 

The study of QCD at finite temperature and finite density is currently 
one of the most attractive topics in particle physics. The heavy-ion col- 
lision experiments aiming to produce the quark-gluon plasma are running 
at BNL and CERN, for which the interesting regime is rather low density. 
C$ ' Moreover a new color superconductor phase is expected in the region of 

low temperature and high density. Numerical study by Monte-Carlo sim- 
ulations of Lattice QCD is a powerful means of investigating aspects of 
the phase transition but the Monte-Carlo method is not applicable directly 
at finite density because the fermion determinant is complex for non-zero 
quark chemical potential /i q . Most of the studies at /i q ^ are done by the 
reweighting method performing simulations at Re(/i q ) = 1,2 ' . 

In this paper, we make a comment on the reweighting method for the 
chemical potential direction. The interesting point is that there are two 
simulation parameters, (3 = 6/g 2 and /i = \i q a. We perform a simula- 
tion at a suitable point (/3o,/xo) and, in order to calculate an observable 
at another point we modify the Boltzmann weight. Then, if the 
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modification factors (reweighting factors) for [3 and fi are correlated dur- 
ing Monte-Carlo steps, there could exist a parameter subspace in which 
the effect of the modification is cancelled and changes of all physical quan- 
tities are small. Therefore it is interesting to investigate the correlation 
between the reweighting factors and to confirm if such a direction exists. 
It is also interesting to discuss the relation between that direction and the 
phase transition line, since we naively expect the physics is similar along the 
phase transition line. Moreover increase of the error due to the reweighting 
could be reduced by the cancellation, which may explain why Fodor and 
Katz 2 could calculate (3 C for rather large /U. 

2. Fluctuation of the reweighting factor 

The reweighting method is based on the following identity: 

i r (Oe AF e AG ) 

Here M is the quark matrix, S g is the gauge action, Nf is the number of 
flavors, F = iV f (lndet M w - In dot M (/lo) ), G = {fi - 0o)P, P = -dS g /df3, 
and Af = f — (/). The expectation value (0}(/3 iAi ) can in principle be 
computed by simulation at (/?o,Mo) by this identity. In this study, we put 
fiQ — 0. If O and e AF e AG — e F e G / (e^ ) fluctuate with a correlation 
during Monte-Carlo steps, (O)m^) has \i and j3 dependence. Otherwise, 
the fi and (3 dependence cannot be obtained, e.g. if e F e G does not fluctuate, 
e AF e AG = 1 and (0}(/3 l([i ) docs not change. Roughly speaking, the differ- 
ence of (0)(^ )M ) from (O)( l g 0; o) increases as the magnitude of fluctuations 
of F and G increases. 

We discuss the correlation between e F and e G . Since e F is complex, 
we separate it into a phase factor and an amplitude. As is shown in Ref. , 
the phase factor and the amplitude can be written by the odd terms and 
the even terms of the Taylor expansion of IndetM, respectively, since the 
odd terms are purely imaginary and the even terms are real at [i = 0. De- 
noting e F = e ie \e F \ and F = £~ i R nfJ , n , \e F \ - cxp{£"=i ReR 2nt , 2n }, 
and e 10 — exp{i^^L 1 Imi?2n-iM 2 " -1 }- We study these correlations in the 
vicinity of the simulation point (/3 , 0). Up to 0[/3 — /3 a , [i 2 ], the reweighting 
factorise^le^le « l + i? 1 ^ + i? 2 ^ 2 /2 + i? 2 /i 2 + P(/3-/3 ). We compute the 
correlations, (A(i? 2 /2)AP), (AR 2 AP), and (A(R 2 1 /2)AR 2 ), which corre- 
spond to the correlations of (e 10 , e G ), (|e F |, e G ) and (e %e , |e F |), respectively. 
Here, (ARiAP) is zero at /x = because R\ is purely imaginary. 
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We use the configurations in Ref. 3 , which are generated by the Nf = 2 
p4-improved staggered action on a 16 3 x 4 lattice. The results are sum- 
marized in Table 1. We find that the correlation between \e F \ and e G is 
very strong in comparison with the other correlations, which means that 
the contribution to an observable can be separated into two independent 
parts: from e* 9 , and from a combination of \e F \ x e G . 



Table 1. Correlations anion; 



d P 



= 16 J 
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(A(Rf/2)AP)JV s T t 1 e 


(AR 2 AP)N- it \ 


(A(Rj>/2)AR 2 )N^ e 


0.1 


3.64 


0.006(29) 


0.312(33) 


0.034(10) 




3.65 


0.059(21) 


0.434(29) 


0.056(10) 




3.66 


0.055(15) 


0.410(26) 


0.022(5) 




3.67 


0.032(15) 


0.397(28) 


0.031(5) 


0.2 


3.75 


0.037(18) 


0.287(26) 


0.029(7) 




3.76 


0.019(10) 


0.353(23) 


0.018(3) 




3.77 


0.037(10) 


0.359(24) 


0.017(3) 



To make the meaning of this result clearer, we consider the following 
partition function, introducing two different fx, fi Q and /i e , 



Z = J J D[/e- Rl ^ + - R3 ^ + -e- R2 ^ +fl4 ^ + -(detM| Al=0 ) Art e" 



(2) 



Then, (ARfAP) = ((ARtfAP) = fg| = 9{xs ~™ s) N site , 2{AR 2 AP) = 

Cpf = ^e^ N sitc at n = 0, where xs and xns arc the singlet and non- 
singlet quark number susceptibilities 4 . The result in Table 1 means that 



gjlgp « gjzgp i i- e - A 1 m the phase factor (fi a ) does not contribute to the 
/3-dependence of Z and \x in the amplitude (// e ) is more important for the 
determination of fi c . Moreover, these correlations have a relation with the 
slope of xs and xns m terms of (3. Since xs — Xns is known to be small 4 , 
this result may not change even for small quark mass. 



Iso-vector chemical potential Next, we discuss the model with Iso- 
vector chemical potential 5 . If we impose a chemical potential with opposite 
sign for u and d quarks: fi u = —fid, the Monte-Carlo method is applicable, 
since the measure is not complex. In Ref. 3 , we discussed the difference of 
the curvature of the phase transition line in this case. Because wc expect at 
T = that pion condensation happens around fi q « m w /2, if we consider 
that the phase transition line runs to that point directly, the curvature of the 
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transition line for Iso- vector /x should be much larger than that for usual /i, 
since m, /2 << ton/3. However, as we discussed above, /i in Eq. 2 does not 
contribute the shift of (i c and the difference from the usual fi is only in 
i.e. \i a — for the Iso- vector case. Therefore the difference of the curvature 
might be small and the naive picture seems to be wrong. In practice, our 
result at small \x using the method in Ref. 3 supports that. Moreover Kogut 
and Sinclair 6 showed that (3 C from chiral condensate measurements is fairly 
insensitive to fi for small fi by direct simulations with Iso-vector \i. 




Figure 1. Left figure is the contour plot of the dispersion around /3 C (0). Bold lines 
show d 2 /3 c /d^ 2 . O = 3.6497, m = 0.1. Right figure shows \/3 c (p) - /3 C (0)\ for real /i and 
imaginary fi at m = 0.1. 



We estimate the fluctuation of the reweighting factor. The phase factor 
e l6 gives a contribution independent from the other parts. As we discussed 
in Ref. 3 , these phase fluctuations cause the sign problem to become more 
severe as \i (or fi Q ) increases. The amplitude of the fermionic part and 
the gauge part are strongly correlated. Therefore the magnitude of the 
contribution to a observable is not simple. We compute the dispersion of 
\e F \e G to estimate the fluctuation. Up to 0\f3 - /3 Q , fi 2 }, ({AO^e )} 2 } « 
fi 4 ((AR 2 ) 2 ) + 2fi 2 (p - /3 )(Ai 1 , 2 AP) + (J3- /3 ) 2 ((AP) 2 ). Then, the line 
of constant dispersion is an ellipse. We write contour line in Fig. 1 (left) . 
Here the susceptibilities of Ri and P and the correlation of i?2 and P 
are computed at the phase transition point, j3 c = 3.6497(16) for m = 0.1. 
We denote the lower and upper bounds of d 2 p c /d/j, 2 = —1.1(4) by bold 
lines 3 . We find that there exists a direction along which the increase of 
the fluctuation is relatively small. That direction is roughly parallel to 
the phase transition line. Because we expect that physics is similar along 
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the transition line, if we consider that |e F |e G is the important part for 
the calculation of (3 C , this result suggests that the phase transition line 
is determined by the quite simple mechanism that the fluctuation of the 
reweighting factor itself is small along the transition line and physics is 
similar on that line. 

Imaginary chemical potential In Fig. 1 (left) , we write also the region 
for fi 2 < 0, i.e. imaginary fi. de Forcrand and Philipscn 7 computed 
d 2 f3 c /dfj, 2 performing simulations with imaginary /i , assuming that (5 C is an 
even function in \i and analyticity in that region. Here, we confirm whether 
the results obtained by real and imaginary /i are consistent or not by the 
method in Ref. 3 . We replace /i by i\x or —ifi and reanalyze for imaginary 
/i. In Ref. 3 , the reweighting factor has been obtained in the form of the 
Taylor expansion in n up to 0(fi 2 ), and the replacement is easy. We de- 
termined f3 c by the chiral susceptibility. The results of |/3 c (m) — /3 C (0) | are 
written in Fig. l(right). Errors are 0(n 4 ). The solid line is the result for 
real \x. The results of [i — > ifj, and [i — > — i[i are dashed and dot-dashed 
lines respectively. The slope at fi — is (<9 2 /3 c /<9^ 2 )/2. We find that these 
results of the slope for real and imaginary /i are consistent. 

3. Conclusions 

We investigated the fluctuation of the reweighting factor. The contribution 
to the //) dependence of a physical quantity can be separated into the 
phase factor and the amplitude. Mainly the amplitude of the reweighting 
factor contributes to the determination of (i c . The fluctuation of the am- 
plitude is small along the phase transition line, which is consistent with 
physics being similar along the transition line, while the fluctuation of the 
phase increases in proportion to fi and causes the sign problem. We also 
confirmed that the second derivative of the phase transition lines at fi = 
determined from imaginary /i is consistent with that from real jit. 
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